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Frequency dissemination and synchronization in free space plays an important role in global navigation satel¬ 
lite system, radio astronomy and synthetic aperture radar. In this paper, we demonstrated a portable radio 
frequency dissemination scheme via free space using microwave antennas. The setup has a good environment 
adaptability and high dissemination stability. The frequency signal was disseminated at different distances 
ranging from 10 to 640 m with a fixed 10 Hz locking bandwidth, and the scaling law of dissemination stability 
on distance and averaging time was discussed. The preliminary extrapolation shows that the dissemination 
stability may reach 1 x 10“^^/s in ground-to-satellite synchronization, which far exceeds all present methods, 
and is worthy for further study. [http://dx.doi.Org/10.1063/l.4921001] 


I. INTRODUCTION 

In recent years, high-precision time and frequency 
signals dissemination via optical fiber has shown rapid 
progress.Dissemination schemes for different topolog¬ 
ical structures, from point-to-point, multi-access, cas¬ 
cade to branching structure have been experimentally 
demonstrated.It is useful for positioning, naviga¬ 
tion, timing, radio astronomy (such as Very Long Base¬ 
line Interferometry and Square Kilometer Array) and 
measurement of fundamental constants. However, the 
fiber-based dissemination schemes can only realize two- 
dimensional frequency synchronization, and the dissem¬ 
inated signals are confined in the fiber network. Conse¬ 
quently, for portable and three-dimensional applications 
or any cases with no fiber sources, free space time and 
frequency transfer is essential and is in great need for im¬ 
provement. Most of the related experiments made use of 
free space laser.For example, F. R. Giorgetta et al. 
developed an optical two-way time and frequency trans¬ 
fer (TWTFT) method which can exchange optical pulse 
between coherent frequency combs, and they achieved a 
residual instability below 1 x 10“^^ at 1000 s.^^ However, 
there are some limitations in optical time and frequency 
dissemination via free space. The transmission quality 
is highly affected by weather conditions, and it is hard 
to maintain an effective transfer link on a large scale be¬ 
cause of the small beam divergence, complex operations 
and environmentally sensitive equipment. As a result, 
optical transfer is not capable of all-weather dissemina¬ 
tion. 

By compensating phase noise actively, we have suc¬ 
cessfully disseminated the L-band frequency signal be- 
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tween RF antennas in free space with a stability of 
3 X 10“^^/s and 4 x 10“^^/day.^^ To study potential ap¬ 
plications of phase compensated radio frequency dissem¬ 
ination scheme, especially the application in the case of 
ground-to-satellite synchronization, we improved the ex¬ 
perimental approach and selected an appropriate field 
test site. The links are set with a fixed 10 Hz locking 
bandwidth to simulate actual ground-to-satellite commu¬ 
nication. The dissemination stability (denoted as a) of 
the links from 0.5 s to 1000 s averaging time (denoted 
as r) at different distances (denoted as d) ranging from 
10 m to 640 m was measured in steady nocturnal environ¬ 
ment. Using a binary regression model, we explored the 
relationship between r, d and a. The regression result 
shows that at a certain averaging time, Ig(cr) degrades 
linearly as lg(d) increases. If this scaling law still applies 
in longer distance, the dissemination stability is expected 
to be better than lxl0“^^/s in a ground-to-satellite syn¬ 
chronization link, which is superior to traditional two- 
way satellite time and frequency transfer (TWSTFT) 
and synchronization methods based on Global Position¬ 
ing System (GPS).^^ 

The details of experimental setup are presented in 
Sec. H, and the experiment results and data analysis are 
presented in Sec. HI. 


II. SYSTEM DESCRIPTION 

In our previous experiment,due to the restriction 
of the experimental site, the local terminal and the re¬ 
mote terminal were located in the same room, and the 
antennas were placed on the roof top of Department of 
Precision Instrument Building at Tsinghua University. 
Consequently, we had to connect each antenna to local 
terminal or remote terminal with RF coaxial cable of 
tens-of-meters length, which introduced undesired phase 
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FIG. 1. (a) Satellite map of the experimental site. The lo¬ 
cal terminal was placed at point O, and the remote terminal 
moved along the testing baseline from point A to point F; A, 
B, C, D, F are 10 m, 20 m, 40 m, 80 m, 320 m and 640 m 
away from point O, respectively, (b) The remote Terminal 
and its antenna under test, (c) The rack mount chassis inside 
the cabinet. 


noise, and restricted longer transfer distance. Recently, 
we built a vehicle-mounted portable experimental system 
to replace fixed equipment. All units of local or remote 
terminal are integrated in a 19-in. rack mount chassis 
with a height of 900 mm, which can be placed in cabinet 
and work with ordinary AC power supply. A Parabolic 
dish antenna with a diameter of 1.5 m, mounted on a 
tribrach with wheels, is connected to local or remote ter¬ 
minal via a 5-m-long coaxial cable. Then the tribrach is 
hitched to a vehicle for mobility. A series of experiments 
were conducted at the Changping site of the National In¬ 
stitute of Metrology (NIM-Changping) on the outskirts 
of Beijing between November 2014 to January 2015. Fig¬ 
ure 1 (a) is the map of the experimental site. Fig. 1 (b) 
shows the cabinet of the remote terminal and the antenna 
under test, and Fig. 1 (c) shows the rack mount chassis 
of the remote terminal, which was placed in the cabinet 
during the experiments. 

The experimental scheme diagram is shown in Fig. 2. 
There are two independent links, which can disseminate 
different frequency signals phase locked to a same ref¬ 
erence clock and recover the reference frequency at the 
remote terminal, respectively. Link / contains Transmit¬ 
ter / and Receiver / while Link II contains Transmitter 
II and Receiver IF Take Link / for example, at local 
terminal the 100 MHz reference clock generates a signal 
Er = Vr cos{uOrtF (pr)- In the slave frequency synthesizer 
there is an oscillator generating Eq = Vq cos(co’t+ 0 o )5 and 
there are several multipliers to generate frequency signals 


coherent to Eq. Here, uj is nominally 100 MHz and (po is 
the initial phase of the oscillator. When the uplink signal 
reaches Receiver /, it is coherent to the output frequency 
signal Ej = Vj cos{ujt po pp) , where pp is a variable 
standing for the phase (with noise) accumulated during 
transmission. Through a coherent transponder in Re¬ 
ceiver /, the signal is returned to local terminal and is 
coherent to Er = Vr cos{ujt (po F 20^), after a round- 
trip. By a two-stage mixing at the transmitter, an error 
signal Ee = Ve cos[(co’^ — uj)t + (0o F (pp — pr)] is gener¬ 
ated, which is fed to a phase-locked loop (PLL), making 
UJ = ujr and F pp = pr F constant when the PLL is 
locked. Hence, Ej is phase locked to the reference signal. 
The same applies to Link II {Ejj). The phase difference 
between Ei and Eu stands for \/2 times the dissemina¬ 
tion stability of each link. For convenience, we ignore the 
factor \/2 in the following discussion. Detailed informa¬ 
tion about the structure of transmitters and receivers is 
described in our previous paper. 

Frequency signals suffer turbulence and temperature 
fluctuation during free space dissemination, which mainly 
result in changes of refractive index. A. N. Kolmogorov 
deduced a structure function of the velocity field in tur¬ 
bulent flows at very large Reynolds numbers: 


DF) = ^0 < 'T < Lo, (1) 

where is the velocity structure constant; Iq is inner 
scale and below which viscous dissipation is negligible; Lq 
is the outer scale over which no kinetic energy is injected. 
This function is often called “2/3 power law.” 

Accordingly, structure function of refractive index is 

Dl{r) = lo^ri^Lo, (2) 

where is the structure constant of refractive index. 
Then after a Fourier transform, we get 

^n(/^) = 0.033C^/^-^^/^ 27r/Lo ^ ^ 27r//o, (3) 

where hc is the spatial frequency, is the power spectral 
density of refractive index. 

The frequency stability of phase compensated round- 
trip dissemination scheme may obey the power law sim¬ 
ilar to D‘^{r). To test this hypothesis and preliminarily 
study the feasibility for ground-to-satellite applications, 
we performed a field test using the 700-m-long baseline in 
NIM-Changping. The local terminal was fixed at point 
O, and we moved the remote terminal to measure the dis¬ 
semination stability (Allan deviation, ADEV) at differ¬ 
ent distances respectively. Note that even perfect phase 
noise cancellation with time delay in round-trip brings 
imperfect phase compensation in one-way transfer. 

A round-trip communication takes 0.15 s between ground 
and medium-earth orbit satellite (MEO) and takes 0.25 s 
between ground and geostationary satellite (GEO). Eor 
these occasions, locking bandwidth of the PLL is as low 
as several hertz.To simulate the ground-to-satellite fre- 
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FIG. 2. Schematic of the portable dissemination system. 


quency dissemination with large time delay and fixed low 
locking bandwidth, we set the PLL locking bandwidth to 
10 Hz. 

The outdoor temperature varied from — 12°C to 
20 °C during the experiments. Low temperature indeed 
brought trouble to the equipment at the beginning. It 
caused some devices working improperly, for example, 
the frequent frequency jumps of frequency synthesizers 
and oscillators. Then we pasted foam heat-insulation 
layer on the inner wall of each cabin in l.dmxlmxlm 
size, making sure that all electronic components work 
properly. However, we can not control the temperature 
in the cabinet, so long-term (more than 500 s) frequency 
stability does not decline as averaging time increases due 
to long-term relatively large temperature fluctuations in 
the cabinet. 


III. MODELING OF SCALING LAW AND DATA 
ANALYSIS 

Figure 3 shows the measured frequency dissemination 
stabilities of the experiment system at different distances 
in steady nocturnal environment. At each distance, we 
measure repeatedly to reduce stochastic error. In the 
shortest distance (10 m) the stability is 2.1 x 10“^^ at 
Is and declines to 5.0 x 10“^^ at 256s; in the longest 
distance (640 m) the stability is 2.1 x 10“^^ at 1 s and 
declines to 2.4 x 10“^^ at 256 s. The frequency stability 
of all the distances is limited to around 1 x 10“^^ after 
1000 s owing to temperature fluctuation in the cabinet. 

The dissemination stability may degrade with distance 
at just like the Kolmogorov scaling law, here 6 is a 
constant. To verify this hypothesis, we establish a binary 
linear regression model, supposing that 

cr(r,d) = 10^T^d^ (4) 

i.e. 

( 5 ) 
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FIG. 3. Measured frequency stability of the dissemination sys¬ 
tem with PLLs closed in steady nocturnal environment. Re¬ 
peated measurement results at each distance are represented 
as the shaded area and the mean values are represented as 
solid lines within shaded areas. 

where z and a are constants. The regression result is 
based on 60 data points, ranging from 0.5 s to 256 s on 
the time scale and from 10 m to 640 m on the distance 
scale, shown in Table 1. Here, SE stands for standard er¬ 
ror. We obtain b = 0.422, a = —0.646, with the 95% 
confidence intervals of b and a being 0.369 ^ 0.475 and 
—0.685 ^ —0.607, respectively. This indicates that the 
dissemination stability declines over averaging time at 
approximate and degrades with distance at ap¬ 
proximate Fig. 4 intuitively illustrates the fitting 

result. The P-value is 0.000, and the adjusted R-Square 
for the whole equation is 0.958, showing great statistical 
significance. Table H displays a strong positive corre¬ 
lation between Ig(cr) and lg(d), and a strong negative 
correlation between Ig(cr) and lg(T). 

In order to explore the relationship between Ig(cr) and 
lg(r) more clearly, we study the scaling law on certain 
averaging times. In this case the regression equation is 


Ig(cr) = a • lg(T) + b ■ \g{d) + z, 
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FIG. 4. Fitted surface plot of distance (d), averaging time (r) 
and dissemination stability (cr) in logarithmic coordinate.The 
black dots stands for the measured data and the plain mesh 
in false-color grids stands for the regression equation. 


TABLE I. Coefficients of parameters in binary linear regres¬ 
sion 



Value 

SE 

95.0% Confidence Interval 

Lower bound 

Upper bound 

z 

-14.072 

0.056 

-14.184 

-13.960 

b 

0.422 

0.026 

0.369 

0.475 

a 

-0.646 

0.019 

-0.685 

-0.607 


simplified into a unary linear equation 

a{d) = (6) 

where 10^ = 10^t“ with a given r. Table III shows the 
regression results at different averaging times. Cl stands 
for confidence interval and R stands for correlation co¬ 
efficient. As r increases from 0.5 s to 256 s, parameter 
b decreases from approximate 1/2 to 1/3, the width of 
6’s 95% confidence interval widens from 0.25 to 0.8, and 
the correlation coefficient also decreases. If temperature 
of the out-of-loop devices and the measurement system 
is well controlled (experimental records indicate signifi¬ 
cant correlation between phase fluctuation and temper¬ 
ature fluctuation with 4^6ps/°C), the difference of 6’s 
regression value between short-term stability and long¬ 
term stability should be significantly narrowed and with 
less apparent drift. 

Suppose the atmosphere is homogeneous with same 
density at sea level, the height is called scale height for 
the atmosphere,which is expressed as 


ksT 
mg ’ 


( 7 ) 


TABLE 11. Correlations of regressors in binary linear regres¬ 
sion 


Pearson Correlation 

Ig(o-) 


ig(T) 

ig{o') 

1.000 

0.425 

-0.883 

ig(d) 

0.425 

1.000 

0.000 

Ig('r) 

-0.883 

0.000 

1.000 


TABLE III. Parameters in unary linear regression 


r / s 

Z 

b 

95.0% Cl of b 

R 

0.5 

-14.004 

0.497 

0.370 


0.625 

0.984 

1 

-14.198 

0.497 

0.319 


0.674 

0.968 

2 

-14.332 

0.465 

0.285 


0.645 

0.963 

4 

-14.508 

0.452 

0.255 


0.649 

0.954 

8 

-14.702 

0.441 

0.230 


0.653 

0.945 

16 

-14.852 

0.410 

0.189 


0.630 

0.932 

32 

-15.030 

0.387 

0.148 


0.625 

0.914 

64 

-15.189 

0.377 

0.078 


0.676 

0.868 

128 

-15.306 

0.365 

-0.027 


0.756 

0.791 

256 

-15.406 

0.333 

-0.065 


0.730 

0.758 


here ks is the Boltzmann constant (1.38 x 10“^^ J/K), 
T is the temperature of atmosphere in kelvins and we 
can take it as 300 K, m is the mean molecular mass 
of air (4.8 x 10“^^ kg), g is the acceleration due to 
gravity (9.8m/s^), thus H = 8800m. The equivalent 
thickness of the atmosphere in the path link between 
ground and satellite varies with the zenith angle, and 
H is the Minimum. For example, the equivalent thick¬ 
ness of the atmosphere from a ground station in Bei¬ 
jing (approximately located at 40°N, 116°E) to a GEO 
satellite at the same longitude is about 13000 m. We 
just take a preliminary extrapolation, substitute param¬ 
eters d=13000m, a=—2/3, 6=1/2, ^=—14 into Eq (4), 
a stability of 1.1 x 10“^^/s, 5.8 x 10“^^/d is estimated. 
The stability of a typical Cesium atomic clock (Microsemi 
5071A) is 5 X 10“^^/s, 3 x 10“^^/d). Consequently, the 
dissemination scheme is capable of transferring frequency 
of a commercial atomic clock without downgrade its sta¬ 
bility. However, as a preliminary study, we have not 
taken into account the path loss of practical ground-to- 
satellite link. The free space path loss from ground to 
a GEO satellite (such as Inmarsat) in L band is about 
190 dB according to Friis transmission equation, and the 
received downlink signal power is about —80 dBm if the 
antenna gain is 30 dB.^^ The signal dissemination and 
processing issues in this case needs further study in the 
future. 
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IV. CONCLUSION 

In summary, we demonstrated a continuous, portable 
radio frequency dissemination scheme via free space. It 
is easy to set up with high mobility, and it can be unat¬ 
tended if automatic control unit is added to the sys¬ 
tem. This flexibility gives it a great potential in three- 
dimensional regional or even global time and frequency 
network and in many scientific and engineering applica¬ 
tions, such as clock-based geodesy and synthetic aperture 
radar (SAR). Besides, a preliminary study indicates that 
the relationship between its dissemination stability and 
transfer distance follows a scaling law, and a stability 
of 1 X 10 ^^/s, 6 X 10 ^^/d is promising for ground-to- 
satellite synchronization. However, long-term frequency 
stability is limited to around 1 x 10“^^ within km-distance 
mainly due to temperature drift in the experimental de¬ 
vices. Consequently, more rigorous experiments focusing 
on longer distance, longer time with fine thermal control 
are worthy and meaningful to verify this possibility. 
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